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AUTOMATED LANDING, ROLLOUT, AND TURNOFF USING 
MLS AND MAGNETIC CABLE SENSORS 


S. Pines 
S. F. Schmidt 
F. Mann 

Analytical Mechanics Associates, Inc. 

SUMMARY 

This report contains a description of the simulation program used to study the 
landing approach, rollout and turn-off of the B73:7-i6o aircraft utilizing MLS and a buried 
magnetic leader cable as navigation aids. The results of the simulation are contained in 
the report and show the concept to be both feasible and practical for commercial type 
aircraft terminal area control. 


INTRODUCTION 


The use of inertial navigations in conjunction with navigation aids such as ILS, 
MLS, TACAN, etc. , in commercial aircraft is rapidly becoming a common state-of- 
the-art technology. As part of the Terminal Configured Vehicle program (TCV), 
sponsored jointly by NASA and FAA, this study was undertaken to: (1) investigate the 
compatability of existing Autoland guidance laws utilizing a Microwave Landir^ System 
(MLS) in place of the conventional ILS navigation aids during landing approach; and (2) 
to extend the concept of automated flight to rollout and turnoff for reduced runway 
occupancy during normal and adverse (CAT III) landing conditions. In order to achieve 
acceptable tracking performance in dry and adverse weather conditions, a buried magne- 
tic runway cable and an aircraft mounted three coil magnetic pickup were investigated. 

For the purpose of the study, a dynamical model of the B737 was generated in 
a computer program (ALERT) to simulate the aircraft forces and moments during land- 
ing approach, rollout and turnoff. A simulation of the existing Autoland guidance 


ix 



equations for the B737 was supplied by the Langley Research Center^ Optimal Kalman 
filtering, utilizing a square root formulation for onboard implementation, and a non- 
optimal fixed gain complementary filter were studied. A rollout and turnoff guidance 
la:w was generated. Both Kalman type filtering and fixed gain complementary filtering 
are shown to be useable with little change to conventional guidance equations presently 
in use in commercial aircraft. 

This report includes a study of the requirements for the Kalman filter and the 
rollout and turnoff guidance laws. 




AIRCRAFr DYNAMICS ^ ^ ^ 

This section contains the derivation; of the equations. used in; the simulation, 

, of the dynamic model of the B-737 during landir^ approach, : capture, decrab, flareij 
touchdown, rollout and tumqff. . .■ t : r; 

A. Choice of the Coordinate Systems 

- ^ i . ? It is customary Vo utilize rinearized 'equatidhS of liiofion desc^^^ ' 
the deviation of the aircraft from a constant-speed, fixfed -flight path angleV wings- ' 
level equilibrium condition in the body-stability axis system. While this is suf- 
ficient for the landing approach phase, it is inadequate for rollout and turnoff due 
to large deviations in air speed, flight path angle, and sizable changes in yaw 
during turnoff. In order to accommodate the broad spectrum of conditions, it was 
decided to use nonlinear equations of motion in a runway Cartesian coordinate 
system to describe the linear acceleration of the aircraft and to use the body axes 
and the Euler angles of roll, pitch, and yaw to describe the angular acceleration 
equations. 

Since the nonlinear variation of the aerodjmamic coefficients is small 
over the range of Mach numbers and aircraft altitudes in landing, and since the 
aerodynamic forces during rollout and turnoff become relatively unimportant when 
compared to the landing gear and thrust forces and moments, it was decided to 
use a simplified form of constant lift, drag, and moment coefficients fitted to the 
B-737 for the landing simulation. Numerical values of these aerodynamic coef- 
ficients were obtained from NASA Langley Research Center and are given in 
Appendix A. 

The orientation and sign convention for the runway inertial Cartesian co- 
ordinate system is given in Figure 1, The aircraft position vector, R(x,y, z) , 
is measured relative to the origin fixed on the center line of the rimway with the 
positive X axis pointing forward along the runway center line. The z axis is 
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Figure 1, Runway (Inertial) Coordinates and Aircraft Position, Velocity Vectors (R, R ). 
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positive down, and the y axis is positive to the right in a right-handed orthogonal 
coordinate system. 


The transformation from the body system to the inertial runway system is 
given in terms of the Euler angles. 

Let 

0 = roll ■ 

0 = pitch 

^ = yaw (measured from magnetic north) 


Then, 


= runway yaw. (measured from magnetic north) 
R 




1 0 0 

0 cos 0 sin 0 

L O' -sin 0 cos 0 


(la) 


where 


T2(0) = 


cos 8 
0 

sin 0 


0 -sin 0 

1 0 

0 cos 0 


(lb) 


Tg(a) = 


cos a sin a 0 

-sin a cos o 0 

0 0 1 


(ic) 


a = ’J' - 


R 


(Id) 


A vector in body system V is transformed into a vector in runway inertial 

system by the matrix L , 

IB 




(2) 
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where 


-hs = 


(2a) 


It follows that the three unit-body-axes vectors a , a , and a maj^ be expressed 

X ^ o 

in the inertial system in terms of the Euler angles (0 , 6 , a) • " See Figs. 2a aind 2b. 
The a^ axis is along the aircraft center line, positive forward. ; 


^ 


cos 6 cos a 
cos 6 sin'icr 
-sin 6 


(3a) 


The a axis is the lateral aircraft axis, poiMing positive along the right wing. 


-sin a cos 0 + cos o sin 0 sin B 1 


^2 ~ I ^ ^ ^ ^ ^ ® ' 

sin 0 cos 9 


(3b) 


The a„ axis is the vertical body axis pointing positive downward; 


sin CT sin 0 + cos a cos 0 sin 0 


^ -cos a sin 0 + sin a cos 0 sin 6, ^ 
dbs 0' cds - 0 ■ 


(3c) 


B. The Aircraft Equations of Motion ■ 

The aircraft inertial velocity vector, R( x , y , z ) is giveri by the tirrie 
rate of change of the position vectdr 


dt ^ ^ 


(4) 


4 






The vector velocity of the aircraft relative to the atmosphere V(v ,v , v ) is 

X ^ o 

given by the aircraft inertial velocity vecfbFand the winds. 


V. = R - [W + G + S} 


( 5 ) 


where 


W 

G 

S 


constant wind vector 

gust vector 

wind shear vector 



The simulation of wind, shear, and gust is not derived in this report and was 
supplied by the NASA Langley Research Center. / . 


The aerodynamic coefficients for lift, drag, and side force sLre obtained 
from wind tunnel data and relate to forces in the so-called aircraft stability axis 
system. In this system, the lift acts normal to the aircraft relative velocity 
vector, V, and in the a , a plane. The drag acts in the a , a plane and 

JL . o X o 

is perpendicular to the lift in the .direction negative to the relative velocity vec- 
tor, V . The side force acts along the a„ axis. To express the stability axes 
(Si* Sg, Sg) in the inertial runway system, we have the transformation 

^IS = ^ 

ij. 

where a is the angle of attack of the aircraft 


or 


a 


. -1 
sin 




V 






A 



V 


V 


(7) 
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Thus the lift acts along the negative s„ axis 

O 


s = -sin a a + cos a a 

o - X o 


(8 a) 


The drag acts along the negative s^ axis 


s j= cos Q! a + sin 0 ! a 

1 1-. O 


(8b) 


The aerodynanalc side force acts along the positive s^ axis 


®2 ^2 


(8c) 


The total forces 'acting on the aircraft, expressed in the runvyay -referenced 
frame, are given by v..- ; 


F = q S L 


IS, 


f-c 1 




f F 1 



XG 


XT 


► + * 

F ■ '' 


F„_ > 



YG 


YT. 

-C " 


'f' 



L 


ZG 


ZT 

^ y 


-(9) 


where 


q = §pV (dynamic pressure) 

S = wing reference area 

Sd ’ ~ effective drag, side force, and lift coefficients 

(see Appendix A) " 

, F_ _ - landing gear forces expressed in the runway, 
. XG , YG ZQj 
■ system 


(9a) 


F , F _ , F = thrust forces expressed in the runway system 
XT Yi -^T 



The time rate of change of the inertial velocity- vector in; the runway 
frame is given by 


-^R = A, 


qs 




dt " ‘*G m ^ S ^ ^ ^ ro i 


YG 


F' ' 

; *ZGJ 


+ — , f cos a a, - sin a a. }; 
m . o 1 o 3 :, 


where 


= < 


0 - 


0 > (gravity vector) 


j*/ . 




k 


magnitude of the thrust 
thrust offset angle insth 
rotation about the ^ axis. 


thrust offset angle insthe. a^:,ao planci Ppsitiye-a jis a- 

X o o 


(io) 


(lOa) 


It is convenient to integrate the; angular. accelerations of the aircraft in 
, ,the body system since "the mass moments of inertia do not change in that frame. 
The kinematic equations , relating the Euler ang^s and thfe body angular rates 
(p,q,r) are given by 






0 = cos 0 q - sin 0 r 

■y"'. 

'Jr = (r cos 0 + q sin 0)/cO8 0 ;.-/K , 5 

0 = ■" p'"+' sin 0 'v ." , ■ 


( 11 ) 


The differential equations for the angular accelerations in the body System 
are glveh'by' ■ ■•y -vy ; ,:v . 


I Coi) + 


( 12 ) 


9 



where 


I 




-I 


XX 

0 

xz 


0 

lyy 

0 


-I 


I 


XZ 

0 

zz 


(12a) 


/I 

w = ^ q ^ 





/ s 

' 

/ s 

T 

- T2(-«) < 

”s 

» + < 

“g 

' M 

^T f 






T ; 


(12b) 



(12c) 



The aerodynamic moments, L_ , M , N , in the stability frame are given in 

O kS O : 

Appendix A. The landing gear moments, L_ , , are in the body frame 

Cj G G , 

and are derived in Section C. The moment due to thrust is given by ^ 


^T 


/ " 
d 


► = < 


T 


T 

i ''t J 


0 

V / 




(12d) 




is the moment arm of the engine thrust and is listed in Appendix A. 


C. Landing Gear Dynamics 

The landing gear is modeled as an airspring and damper. Strut com-^ 

pression is initiated whenever the strut load exceeds a preset preload, P . 

Li 



For loads less than the preload, tire compression provides a linear spring (initial 
tire compression). Elastic deflection of the strut and inertial acceleration of the 
unsprung mass of the strut and tire are considered high frequency phenomena and 
are neglected in the analysis. 

In order to determine the forces and moments acting on the aircraft due 
to ground reactions in landing and braking, the locations of the main gear and 
nose strut with respect to the center of gravity are required. The locations of 
the three extended, uncompressed gears of the B-737 are given in Ref. 1 and 
illustrated in Fig. 3. The positions of the three gears are measured from an 
origin at the leading edge of the MAC. 

*th 

Let the coordinates of the i gear with respect to the leading edge of 
the MAC be given by 

' ' X. , y. , ¥. (see Table i) ' 

i -^i i ' ' 

i ■ ’ ’ . I • ' . 

During landing, rollout, and turnoff, it is assumed that the aircraft experiences 
small values of pitch and roll. Large values of yaw are permitted. Thus, small 
angle theory is used to model roll and pitch motions but the full nonlinear theory 
holds in yaw. 

The height of the gear above the runway . is given by 
liiG = " ^iG " ~ ^ ~ ®^^i ^ 

where 

% 

g =? center of gravity expressed- as the percentage of mean 
aerodynamic chord 
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(30.5") 

. 7747 m 



CG AT 
LEADING 

edge mac 


LANDING GEAR LOCATION 


CG 



EXTENDED GEAR HEIGHTS 


Figure 3. Landing Gear Geometry 



TABLE 1 


EXTENDED GEAR LOCATIONS 


(Distances in meters from c.g. ) 



X 

y 

z 

Gear No. 

Right Main Gear 

-1. 8288 

2.6152 

2.9809 

1 

Left Main Gear 

-1.8288 ; 

-2.6152 

2.9809 

2 

Nose Gear 

8.6258 

0 

2.6731 

3 



Ground contact occurs whenever 


h <0 or > 0 

iG iG 


(13a) 


Since the strut is preloaded, no strut deflection occurs until z exceeds the 

ICj 

preload tire deflection, . Thus, for 


0 ^ z s z. 
iG ipL 


(13b) 


the strut load is a linear function of z.^. 

iG 


F. 

1 



'iG 


ipL 


(14) 


Once the preload has been exceeded, the strut begins to compress and the force 
is given by the isothermal compression law 


^iG ^ipL ] 

1 . 

i 

where f,. is the bottom-out maximum deflection permitted by the strut design. 

Equation (15) holds for all z such that 

iG 



(15) 


z, ^ 

ipL 


iG i 


(15a) 


The strut deflection is given by 


z_ = z.„ - z. ^ 

iD iG ipL 


(15b) 
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The damping force is generated by the strut deflection rate, , 
given by 

ZjQ = z + Py. - q(x. + eg*) 

The damping force of the strut is given by 


^iD “^iDhiol'^iv 

where c.^ is the damping coefficient, which is modeled as a function of the 
strut deflection 


c. 

IV 


a z + a z 
li iD 2i iD 


where a , a are strut damping design constants. Figures 4 and 5 illustrate 

•L I ^ I 

the strut compression curves for the main and nose gears. 


The net vertical strut force in the body frame is given by 


F 

Ni 


F. + F_ . 
i Di 


The tire drag force may now be computed as a function of the normal 
strut force and the friction coefficient. The drag force acts in the plane of the 
tire and in a direction opposite to the axle velocity. We have 



The tire friction coefficient, fj,, is the sum of the rolling friction coefficient, ^ 
and the friction coefficient developed by brake application, . We have, for 
dry runway conditions. 


( 16 ) 


(17) 


(17a) 


(18) 


(19) 


15 





static Strut Force (Newtons) 





.5 


( 20 ) 


pmax 


where b is the brake pressure and b is the maximum pressure available, 

p pmax 


For wet runways, we have (from Ref. 1) 


+ 100 

^ G 

10/ Vr +;200 b 


pmax 


where is the ground speed in knots. The tire friction curves are given'in^ 

G ./ 

Fig. 6. ‘ i 


The guidance law for brake pressure application is given in the Guidance 


and Control Section of this report. 


Finally, we have 


M I 

' i 

wherje is the coefficient of rolling friction. Braking is not available. onTlie . 

nose wheel so that only rolling friction prevails. 

The landing gear lateral forces and steering is taken from the theory 
developed by W. B. Horne in Ref. 2. We first compute the tire compression 


coefficient v . : 


I^Ni 


'i “ ^ 

( p. + . 08 p . ) w. ( w. d. ) ® 
i ri 1 i i 


T7. > . 09 


18 




8 

:y (kts) 


9 







the tire compressipn, is given by 


6.. = w.(.03 + .42 Tj.) 

ti i' 'i' 


(22a) 


If 


Tj. < . 09 


tire compression is - - n 


6.. = 7 W. T7. 

ti 4 1 'i 




(22b) 


The tire cornering force coefficient is given by 




(23) 




( ^hPj factor, 2 •, .multiplying ir /l Sp 

•'.i-.* ■. i. 

Thus.^ for, a S;ingle tire strut, the . factor would^b I,.).., 


For 






fi.; s . 0875.d, 
tl •' I 


N. = C... 

i . Ni 


^ ‘ 6 

. 067 + . 34 ^ 
a. 

' I 


(24a) 


and for 


6,. < . 0875 d. 
tl i 


N. = C^T- 
i Ni 


6 ,. 

1 2 — 

* d. 

i 


8 . 8 




(24b) 



The lateral tire cornering force acts normal to the plane of the' tire 
and is given by 


Pr • ) 

Lt 11 I 


where is the ground track angle 


£. = tan 
i 


y + r C( X. + c g )cos ct - y. sin cr ] 


(cr = 'I' - ) 


and V. is the steering angle of the wheel. For the main gear, ■>: = 0 ; and 
for the nose wheel, y. is supplied by the control logic. 


The plane of the main tire is paraUel to the a , a body frame. How- 

' , X ■ O' 

ever, the nose wheel develops a drag and cornering force which must be re- 
sblved into the body axes/’ Let y ■ , F ^ , ahd F_ J be.the nose wheel steering iv 
angle and horizontal forceis developed in the plane of the nose wheel arid perpen- 
dicular to it, Then we have, for the forces in the body system. 


F = F - 0F " - 6 ' y ’F 
Xi Hi . Ni °i3 ^3 L3 


F = F + + fi yF 

Yi Li Ni n3?3 H3 


F • = 6F - 0 F 4 F - 6 y (6F + 0F ) 

Zi Hi Li Ni ®i3^3^ L3 ^ H3 ' 


To obtain the moments in the body frame, we have 


21 



The lateral tire' cornering force acts normal to the plane of the tire 
and is given by 


Fy • ^ N.(y. + a - e. ) 

Li 11 i 


( 25 ) 


where e. is the ground track angle 


e. = tan 

i 


^ y + r ^ X. + c g )cos a - y. sin a ] 


; ‘ : (25a) 


(a = -if - 


and Vj is the steering angle of the wheel. For the main gear, y. = 0', and 
for the nose wheel, y. . is supplied by the control logic. 


The plane of the main tire is parallel to the a , a body frame. How- 

J. O' 

ever, the nose wheel develops a drag and cornering force which must be re- 
sblved into the body axes ^ Let y„ , F„ , and F_ ■ be the riose wheel steering 
angle and horizontal forces developed in the plane of the nose wheel arid perpen- 
dicular to it. Then we have, for the forces in the body system. 


F 

Xi 

F - 0F - 6 y' F 
Hi Ni °i3 ^3 L3 



F + <t)F + 6 V F 
Li Ni °i3^3 H3 

(26) 

'■zi ' 

6F - 0F + F - 6 y (0F + 0F ) 

Hi ^ Li Ni ®i3^3' L3 ^H3 ^ 



To obtain the moments in the body frame, we have 


21 




Z ^Yi^iG^ 


i=l 


IVL 




i=l 


(27) 


"b = 


I 

i=l 


[Ey.(X| + cg )- 


To obtain the forces in the ninway frame, we have 


XG 


" ^3 ^L3 ^ Z ^Hi 

cos cr - 

^3 ^H3 ^ Z ^Li ] 

i=l 


i=l > 


[ ^H3 Z 

Ft • 
Li 

sin a + 

" ^3 ^L3 Z ^Hi] 

‘ i=l 



i=l / 


o 

^ZG " Z ^Ni 

i=l 


vci 
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Figure 7. Navigation System Block Diagram 








where is the specific force measured in the body frame by the acceler- 

ometers. Since the accelerometer and gyro data are sampled at 20 Hz, we 
may integrate Eq. (30) with a simple Taylor series 

R(t + At) = R(t) + R(t) At + R(t) ^ 

• • •• 

R(t + At) = R(t) + R(t)At 

At = . 05 sec. 

In the presence of error sources cited earlier, the dead reckoning esti- 
mates become useless after some time, and navigation aids and filter techniques 
are required to correct for the errors. In the subsequent sections, we describe 
two such filter systems. The first is the Kalman filter. This filter utilizes a , 
priori estimates of the uncertainty, in the error states to be corrected and the 
expected variance in the observations, to obtain a minimum variance estimate 
of the correction to the states as a linear function of the difference between the 
observations and the expected value of the observation. 

The second filter system is called a complementary filter. It uses fixed 
gains to estimate the corrections to the position and velocity of the aircraft as a 
linear function of the observation residual. Since it knows nothing about the a ' 
priori uncertainty of the aircraft errors or the observation noise, it depends 
solely upon the concept of exponentially correlated random noise to smooth out 
noise. As shown in the results of the study, the coniplementary filter does as 
well as the Kalman filter in landing approach, but does not do as well in turnoff 
because of the effect of gyro misalignment biases. 



B. Kalman Filter Implementation : , . 

, The optimal filter algprithm involves a very large number of machine opera- 
tions. The basic algorithm separates into the following: , ; ' 

1. Calculations associated with a change in the time reference 

of the filter, and •' 

2. Calculations of the incremental state estimate from the 

measurement information, ' 

In order to. keep the number of opera.tions within the capabilities of current air-r 
borne coinputers,' the accelerometers, and gyros are, sampled, at a rate of . 20 Hz , 
the observations are sampled at a rate of 10 Hz , and the filter is, updated at a 
rate of 1 Hz , 

Past experience has shown that the square root implementation of the opti- 
. , ; mal filter algorithm (Ref. 3) can reduce the effects of numerical errors to insig- 
nificant levels. The square root implementation is therefore incorporated in the 
proposed design. Modeling errors are compensated for by the appropriate use 
of random forcing functions, , This technique causes the more recent measurements 
to be weighted higher than past measurements;, therefore, the estimate tends to 
follow the more recent measurements. , 

The proposed filter will operate in a mariner illustrated in the sketch below: 



At the start of the sequence, we assume that the filter has its covariance 

matrix referenced to time t, . At the times t, , t, , , t„ , measurements 

k k 1 9 
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are accepted and residual sums and partials are computed and saved in pre- 
processing routines. After accepting the measurement at time t , the residual 
sums and partials are transferred to arrays for processing by the filter j and the 
locations used for the pre-processing are cleared for use in pre-processing the 
tj^^^ and subsequent measurements. The residual sums are processed by the ’ 
filter and the incremental state change computed in lower priority logic. When 
these lower priority calculations have been completed, the logic sets markers 
which cause the higher priority logic to pick up the state change and add it to the 
system. Meanwhile, the lower priority logic updates the covariance matrix to 
the time t^^^ and readies the filter for processing the residual sums in the next 
interval. The onboard program operations for executing such Ibgic will be de- ■ 
scribed in a later report, - : ? iy 

The error state is assumed to obey a linear differential equation 

dx'=- F' dx + F' ri ■ :i;ra (32) 

X n ' 

\ 

■ dx = the n-component error state vector (n = 15) ’ 

F ' = an nxn matrix v j o: 

X 

F = an nxm matrix ' ’ 

n 

n = an m-vector of random forcing functions for 

compensation of error growth caused by unmodeled 
error sources. 

An approximate solution of Eq. (32) is 


where 







(33) 
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where 


$ = the transition matrix 

= the forcing function sensitivity matrix 

u(tj^) = a constant ( in the interval tj^ to vector 

for approximating the effects of the noise vector, 
n , of Eq. (32) 


' ' The covariance matrix for the error state at time tj^ is giveh by 

^ ' ^ ■ P(tj^) - W(tj^) wctj^)"^ - e ■ (34) 


where 


W(tj^) 


the Square root covariance. ( is balculatM in 
the proposed square root implementation of the 
optimal filter) 


£ [ ] = the expected value operator. 


v^i7. 


We assume that u(tj^) is a random independent vector such that 




i 7^ -t 


T . 


where U(tj^) is the square root of the process noise variance. 

The appropriate use of the expected value operator with Eq. (33) gives 
the time update of the covariance matrix 


(35) 




T T 


xF 


( 36 ) 


,u 




T 

The matrix, of Eq. (37), has dimension (n+m)xn. The 

Householder algorithm described in Ref., 4 can be used to reduce this matrix to 

an upper triangular form; that is, all the terms below the diagonal are zero in 

T 

the reduced matrix. The matrix reduction algorithm leaves the product W W 
invariant. 

The measurement residual. Ay , , is defined by , 

Ay(t) = y(t) - y(x, t) , . , (38) 

where y(t) = measured value 

y(X;,t) T= the computed va^e of the measurement based, , , . 

on the current estimate of state indicated by 
the output of the navigation equations. 


It is assumed that the residual is related to the error state by 

y(t) = H dx(t) + q ' (39) 

where 

q = the random error in the measurement 

Since we may have an estimate of dx given by dx from other measure- 
ments, we write Eq. (39) as 
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(40) 


Ay = y(t)-Hdx(t) = Hdx(t)+q 

= our best estimate of the residual including 
any past measurements whose effect is not 
yet included in the estimate of state carried 
by the navigation equations 

As was mentioned in the beginning of this section, the filter operates with 

the reference time t^^ for measurements in the interval t^^ to t^^^^ . This 

means that we estimate dx(t ) rather than dx(t). Hence we make the assump- 

tion for t, < t < t, , that 
k k+1 

dx(t) = $(t;t|^) dx(tj^) (41) 

Equation (40) can then be written as 

Ay = y(t) -H 4(t;tj^) dx(tj^) (42) 

One may note that Eqs. (41) and (33) are in disagreement since the in- 
fluence of the random forcing function, . is ignored in Eq. (41). The 

effect of this simplification is that measurements are not quite optimally weighted. 
Experience has shown that this simplification is justified if the batch interval 

A. = t, , - t, is short compared to natural periods of INS error growth (e. g. , the 
k k+1 k 

Schuler period). 

We let 

= W(t^) Wi/ hJ /(H^ W(t^) +Q) (43) 

Q = e[q^] 
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Then the estimate of the error state following inclusion of the measurement is 

d x( t ) = d x( t ) + K Ay (44) 

k a k b m ' ' 

In Eq. (44), the subscript notation ( ) and ( ) means before and after inclusion 

D SL 

of the measurement, respectively. 

T 

The square root covariance matrix, W(tj^) , after inclusion of the 
measurement is given by 

W(t^)J = W(tj^)^- (45) 

where 

T T 

^ «m 

V = 

X = ( c + Q) (i + /qAc'^C + Q) ) 

- Equation (45) is referred to as Potter's algorithm (see Ref. 3). 

In order to minimize the number of measurement operations, we will use 
accumulated residuals over about a one-second period. This means that the 
residual, Ay(t) of Eq. (38), is the accumulated residual for all the measurements 
of the same type in a half-second interval. The partials of the observations with 
respect to the error state, H of Eq. (43), are calculated and accumulated simul- 
taneously with the residual accumulation. The Q of Eq. (43) is the assumed 
, variance of the random noise in the accumulated residual sum. 
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C. Error Model of the Estimated Variables 

We choose for our aircraft error state a vector of 15 variables. These 

are: 

R error in position (3) 

R error in velocity (3) 

U) gyro drift rate (3) 

y bias error in observables (3) 

z vertical accelerometer scale factor error (1) 
w error in the horizontal wind (2) 

These prove sufficient to maintain a stable continuous estimate of the aircraft 
states and to provide the Autoland guidance equations and pilot displays with 
usable information over long time periods (several hours). 

The time history of the 15-dimensioned error vector is described by dif- 
ferential equations. The derivative of the position error vector' R(t) is given by 

^ R(t) = R (46) 

Let the specific force in the runway frame be { f }j, which i si given by 

(ni= e. ' (47) 




( 48 ) 


\ 



The differential equation for the error velocity vector . is given by 

- ; V ‘ ^.5. =. [f']j+:Xe3g, - (49) 

where { £ ]* is white noise. 

R 

The three gyro drifts, the three sensor biases, the horizontal wind com- 
ponents, and the vertical accelerometer scale factor are described as exponentially 
correlated random variables, and their differential eqiia:tibns are- 


dt 


X. 


i 


— X. + a 
T. I 
t 


,/r 

I V T. 


i = 7,15 


(50) 


where r. is the correlation time constant, a. is the variance, and e. is 
white noise. 


. ( . The, soluU of the above ec^uations, for the short time span of the filter 
update, is given by 




At 


R(t + At) = R(t) + R(t) At + { f 


.. R(t + At) - R(t) + { f LAt + { £ }~ At 

‘ ‘ .'1 X. -/v' Xv v;;. ■: ‘ 'i 

tJ.(t + At)) = {e„ ‘ Sj(t)} >,{<7, 


'i 


(51) 


1 = 7, 15 

Equation (51) provides the elements for the state transition matrices, 

and ♦ (t. , ,t, ) . These are listed in Tables 2 and 3. 
u' k+1 k' 
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TABLE 2 


STATE TRANSITION MATRIX 
Nonzero Elements 


* =1 

1,1 




„= f 


At 


* = - f 

^■1 /\ 


At 


1,8 z 2 ^1,9 y 2 


* =1 

2,2 


$ =1 
3,3 


* 2 , 5 - 


* 2 . 6 = 


4 = -f ^ 

2,7 z 2 

.2 


= f 


At 


3,7 y 2 


& = f 

* n rv 


At 


2,9 X 2 


4 =.f 4, 

3,8 X 2 3,15 2 


A= ^ 

4,4 


«>. = f At 

4, 8 . z 


„= -f At 

4,9 y 


$ =1 
5,5 




*5.9= 


$ =1 
6,6 


4»„ „= f At 
6,7 y 


*6.8= -'x^‘ 


$ = At 

6,15 


-(At/Tj -(At/rj -<At/rJ 

& . = e $=e $=e 

^ 7,7 ^ 8,8 ^ 9,9 


®.;=e $ =e $ = e 

10,10 ^ 11,11 12,12 


-(At/r ) -(At/r ) -(At/r ) 

$ =e $ = e ^ $ =e 

^13,13 ^14,14 ^15,15 
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TABLE 3 


ROOT PROCESS NOISE MATRIX 
Nonzero Elements 






D. Kalman Filter Data Updates 


The guidance and control functions require smooth and current data cor- 
rected for the errors listed above. Since the Kalman filter updates the aircraft 

corrections every second, AT , and the dead reckoning estimates of Eq. (31) 

s 

are updated every . 05 seconds, AT, the corrections to the dead reckoning 
estimates are smoothed and updated as follows: 


Let the Kalman correction vector of R(t) and R(t) be Ax(t). We de- 
fine an exponential smoother, a six-element vector function, CS(t), to be given 
by 


-(AT/AT ) 

CS(t + AT) = e ® [CS(t)- Ax(T )] + Ax(T ) 

s s 

with CS(0) = 0 


(52) 


The corrected, smoothed values of the R(t) and R(t) vectors are given by 


R(t) = R(t) + CSj^(t) 
R(t) = R(t) + CS^(t) 


(53) 


To correct for the gyro drifts, we have the following: 

A ■ A. , ■ 

Let Ap , Aq , and Ar be the current best estimates of the gyro drift 
rates and let 0, and 0 be current readout of the uncorrected Euler angles. 
Then, to obtain the best estimates of the Euler angles, we have 
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(54) 


0(t) = 0(t) + [COS ^(t) Ap(t) + sin A"q(t)]/cos 6(t) 

6(t) = 6(t) + cos ^(t) Aq(t) - sin '^{t) Ap(t) 

4^(t) + Ar + [0(t)- 0(t)] sin 0(t) 

d(t) = 

Finally, to obtain the best estimate of acceleration vector in the runway 
frame, we have 


• • 

0 


f ' 

0 

R(t) = < 

0 

G 


0 

. ^2, 


where Az is the best estimate of the bias in the vertical accelerometer. 

E. Complementary Filter 

A complementary filter can be designed to provide a set of nonvarying 
gains, without any relationship to the uncertainty in the aircraft state or the 
biases in the INS and MLS measurements. By adjusting the complementary 
filter gains to pass the relative low frequency of the true aircraft acceleration, 
it is possible to filter out the high frequency noise in the MLS measurements. 

The main criterion for such a filter is to provide a stable estimator, one that 
decays disturbances exponentially and lets the low frequency motions persist. 

A block diagram of the complementary filter is shown in Fig. 8, 

Let one of the coordinates of the R(t) vector be x.(t) . The differential 
equations of the complementary filter for the x. coordinate are given by 
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Figure 8. Compiementary Filter for Navigation With Body-Mounted Accelerometers Using MLS Measurements 








( 56 ) 


dt 






d . 
dt i 




where 


— til ~ 

x.(t) = i component of R(t) [ Eq, (30)] 


It Ls plain that if the residual is zero, the complementary 

filter will provide the dead reckoning solution of the aircraft motion in the runway 
system. Moreover, the system is linear with constant coefficients K , K , and 
K (the filter%ains), ‘ with nonhomogeneous forcing functions due to the noisy MLS 

O']-.. 

measurements, and the aircraft measurement of acceleration obtained from the gyros 
and accelerometers. .By choosing the values of , Kg , and Kg , the solution 
can be designed to be stable and to filter out the high frequency noise. 


; Since the MLS measurements are available only in a discrete form at 

fixed intervals and the output of the accelerometers are available at fixed inter- 
vals, we may, produce a pseudo MLS measurement which is continuous and is 
equal to thei actual MLS measurement at the discrete sampling time, 


We .model the pseudo MLS, to be: 

- - 


We choose x,,.. „ (t ) so that at t = T , the observation time, the pseudo MLS 
MLS.' o' 

is the same as the observed MLS. 
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=SlLS <‘o> ■ •'5*. (W - *.,(‘)AT - [5,(T) +x,,(tjl 
i 


obs 


i2' o' 


AT^V-"^'" 




With this, definition of the MLS measurement, ^we can effect, an explicit solution 
of the system of differential equations [ Eqs. (56) ] which has the required con- 
dition of giving the dead reckoning solution if the residual is zeroi and which 
corrects the state to match the MLS observa.tions without the high frequency 
component. The solution is given by 


AT* 




*12™ = *12< *0 > " " *13< ‘o ' 3 " ¥ =‘mLS.“o ) 3 V = 

1 


*13™ “ *13‘*o>^‘'3[’‘MLS.'y -^l''o>^ 

I 


where the constants b , b , and b„ are functions of K , K , K_ , and AT . 

X ^ u , 1^3 


1 


a + [w^+ 2a^]e“^^'^ cos(u>AT) 


v.’+\ 


0!^- j3(/S^+w-j ]e sin(U) AT)/(*J 


2 . 2 . 




2a^fl[e e ^^-cps(wAT)] \ 1 

+ [ W^) - (U)^+ p^)' -] e‘^ sin( CO A T)’/w | a + (/+ 


"3 = ^3 


a[e COS(WAT)] \ 


- [co^+p^-a^le sin(coAT)/to 


2 2 2 
a +{w - 2a fi 


(57a) 




... (58) 


■^-2 a /5 


(59) 

2ap 



The constents, a, p, co are the roots of the characteristic equation of the 
system defined by Eqs. (56) 




( S + /!!)( S + ^-ia;)(S+/? + iti^) = S^ + K^S^ + KgS +Kg = 0 


= a +2p 

2 2 

K_ = 2aP + e +W 


(60) 


Kg = Q!(fl^ + a>^) 


By choosing small positive values of a, P, and co we can obtain a stable low 
pass fitter.; We recommend values of . ' • ^ 


a = .08 


(61) 


ot 




,, B = tu =, 

VF 


To provide the best estimates of the aircraft position, velocity, and ac- 

S'. , S:' V- 4 . , ' ' 

■ delegation in the runway frame using the complementary filter output, we have 


Xii(T) 

R(t) = <|x (T)^ 


R(t) = < 


Xi2(l) 

"32<^> 


(62) 


Xig(T) 

R(t) = ^ Xgg(T) ^ +R(t) 
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It is to be noted that the complementary filter provides no cbrrection to 
the errors in the estimates of aircraft attitude due to gyro drift or to the position 
error due to MLS biases. •, s 





m. 


'THE MLS OBSERVABLES 


The characteristics of the MLS observables are described in documents 
supplied by NASA Langley Research Center • The system used is for 
conical coordinates of azimuth, elevation, a, and range. The range and 
azimuth origin is at the center of azimuth antenna facing outward toward the ap- 
proaching aircraft along the runway centerline. The elevation origin is offset 
from the runway and is at the center of the elevation antenna. Thus, the Car- 
tesian coordinate system for the MLS observable is the negative of the runway 
inertial system in the x and z directions and in the same direction as the y 
of the nmway. We have, in the rimway coordinate system 


X ‘ 


x' 


X ’ 


y 

z 

MLS 

-y 

z 

+ 

Runway 

-y 

z 

(63) 

Antenna 


where the antenna vector specifies the azimuth antenna position in the runway system. 
The origin of the MLS system is at the azimuth antenna and the position of the 
elevation antenna is measured relative to it. (See Fig. 9a) 

Let the MLS coordinates of the aircraft be X,, , Y,, , Z,, . Then the 

M M M 

observables in the MLS coordinate system are 

RANGE = JX^J + yJ + zJ 
V M M M 

/3 = sin“\-Y^ ^RANGE) (64) 

_1 ^M " ^OE 

= sin — ■ — — _ . ■ - — — 


a 



AmC3RAFT 



GROUND' I 

PRGJEGTION,^ j 


Figure 9(a). The MLS Coordinate System. 
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In the above, X and Y are the coordinates of the elevation antenna from 

Uiii 

the azimuth origin, and Z is the coordinate of the elevation antenna in the 

OE 

MLS system. 

To invert the solution and obtain X,,. Y,, , Z,, from RANGE, - 6 , and 

M M M , 

O', we have 


= g +/g^-h 


Y^ = - RANGE ( sin ) 


(65) 


wliere', 


M 


= y.(RANGE)2-xJ-Y^2 


" ■ ® "" ^OE ^ ■ 

h = (X^g + Y^^^)(sina)^ + yJmRANGE)^ ( cos a)2 - 2Y^ YQ^ sin af 

It should be noted that the above equations do not account for height differences of the 
elevation antenna from the MLS origin which has been defined as the coordinate 

OE 

Finally, to obtain the aircraft coordinates in the runway inertial system, we have 


X 


X 







M 

y 

— 

y 

•- 

”^M- 

z 

Runway 

z 

Azimuth Ant, 

- ^M'- 


( 66 ) 
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A. 


MLS Noise Model 


The noise characteristics are obtained from Ref. 5. For the exponential 
correlated noise of RANGE, |S, and a, we have 

"r<*> ■ “(‘)^Ag>/Az(t-l) 

where u(t) is a random number of zero mean, a. is the standard deviation 
of the i^^ measurement and 



-(AT/T.) 
A. = e ^ 

i 

(68a) 

For initialization 




rj.(O) = u(0) 

(68b) 


Table 4 contains the values of the noise characteristics used in the MLS error 
model. 

To simulate data dropout, we use a random number with uniform distri- 
bution from -1 to 1 and set the observation equal to zero whenever the random 
number exceeds . 98 . Thus, 2 % of each observation string is lost. Since there 
are three observation types, 6% of the pseudo state observations are lost. This 
percentage can be varied. 



TABLE 4 

MLS ERROR MODELS 



Note: IGRS = Independent Gaussian Random Sample 

Initialization is achieved by setting y(0) = u(0). 

See Figures 9(c) and 9(d) for illustrations of MLS jitter in azimuth and 
elevation. 






To simulate bad data, a similar random number with a uniform distribu- 
tion is used, and an average percentage of the observations is set equal to a 
large number ( 1000 times the. bias) Wen though the data valid flag is set true". 
Thus bad data enters’ the system marked good a small percentage of the time. 

Biases are included in the simulation. 'Values of the biases are initially 
chosen from a random number generator and retained for the simulation run. 

The biaS: values are also given in Table 4. 

MLS "jitter" or "stagger" is simulated to produce the data train sequence 

as shown on Figs, 9b, 9c, and 9d. 

' ■ ' . ■/ ' 

To accommodate data dropout, /MLS invalids, and bad data generally, 
the recommended course of action is to replace the bad data type by the air- 
craft estimate of the data missing (or considered rejectable) by the best esti- ' 
mate of the data type based on thfe: current values of the vehicle state. Thus, if 
data is missing, we replace the bad data by range, azimuth, or elevation from 

A . A A 

Eq. (64), using , and ' ~ computed from Eq. (63). This mode is 

•i 

still under study. 

. . ■ I' . ■ ' 

B. Radar Altimeter I 

The radar altimeter readout is activated when the aircraft estimate of 
altitude is below 150 ft. The MLS elevation readout is simultaneously de- 
activated. For the purposes of this study, the error model of the radar al- 
timeter is described by an exponentially correlated bias variable driven by 
white noise. Thus, the radar signal readout 4s given by 
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SU&-SEQUENCE#V: 



0 10.2 


26 31 36.2 38.6 53.8 53.0 64.0 ms 



WORD 

#1 


NOTES: (1) AUXILIARY DATA (1 WORD) OR BACK ELEVATION 
tt) 360® AZIMUTH 


Figure 9(b). MLS Conical Scan Time Sequence. 
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where 


h(t) = - z(t) + bj^(t) 




(69) 


(69a) 


C. MLS Partials 



The Kalman filter requires the partial derivatives of the MLS observables 
■with respect to the 15 estimated variables. - 


We have, for the azimuth observable. 
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For the bias in azimuth, we have 


b/S 


= 1 


bb. 


R . 


(70a) 


For the remaining 11 variables, we have 


bx. 


= 0 


(70b) 


The elevation observable is 


bg 

bx 




bg 

by 




^oe" 




(71) 


bg 

bz 






-1 + 




R 


OE 


where 


^OE ^OE^, "^OE^ ^ ^OE^ 


(71a) 


For the bias in elevations, we have 


bg 

bb 

g 


(71b) 
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For the remaining 11 variables, we have 


ba 

bx. 


= 0 


The range observable is 


b RANGE 
5x 


X 


M 


RANGE 


6 RANGE 

by 


\ 


RANGE 


' bRANGE 
bz 


M 


RANGE 


For the bias in range, we have 


bRANGE 

bb 


= 1 


For the remaining 11 variables, we have. 


bRANGE 


bx. 


For the radar altimeter, we have 


bh(t) ^ 
bz 


For the altimeter bias, we have 


bh(t) 

bb. 


= 1 


(71c) 


(72) 


(72a) 


(72b) 


(73) 


(73a) 
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For the remaining 13 variables, we have 



IV. 


MAGNETIC LEADER CABLE SENSOR 


While the MLS sensor described in the previous section proves adequate 
for landing approach through touchdown, problems in multipath distortion and 
other biases create serious problems against relying solely on the use of MLS 
in rollout and, turnoff . In prder to provide a navigation sensor aid with smaller 
errors, the concept of a single magnetic leader cable embedded along the center- 
line of the runway and along the turnoff lane is recommended. 

Initial studies of such devices began. in the 1960's. Reference'6 contains 
a description of a magnetic cable system and a three-coil sensor developed in 
England. A series of studies for use in automotive steering along reinforced 
concrete highway systems is presently sponsored by DOT. Reference 7 contains 
a description of such an automotive steering aid for lateral control. 

,1 : 

For our purposes, the British system and experience seems to be more 
applicable and is incorporated into this report. 


A. Three-Coil Sensor 


The magnetic field is generated by a single conductor set beneath the 
nonconducting concrete and along the runway centerline. A detector, consisting 
of three coils having an equal number of windings of conducting wire with iden- 
tical characteristics, is used to generate three voltages. The coils are mounted 
about a spherical, nonconductive shell in three perpendicular planes. Figure 10 
illustrates the leader cable and the three coils. 





Let the Cartesian components of the magnetic field be 


H = (H , H . H ) 
X y z' 


( 74 ) 


as illustrated in Fig. 11. We have, for the three induced voltages. 


V, 2: H 
. 1 z 


v„ s; H cos ( A'l') 

2 y 


(75) 


V s H sin (A'i') 

•j y 


For the yaw deviation, we have 


tan 


(76a) 


and for the lateral displacement 

H 

Ad = Az 

rl 

y 


Sf Az 


y 


2 2 
^^2 


(76b) 


The error in yaw and the error in lateral displacement are directly ob- 
tainable from the voltages. To obtain the time rates of these quantities, which 
are required for damping stability in the Autoland guidance equations, we have 
chosen not to rely on differentiation of the magnetic coil signals, but to generate 
the required information analytically. The equations for these functions are 
given in the Guidance Section of this report. 




B. Error Model of the Magnetic Coil Sensor . 

-.fr •, 

The distortion in signal strength is due to several sources, as shown in 
Refs. 6 and 7. 


1. ground conduction > ' ; 

2. nonlinearity 

3. phase modulation ’ ■ ^ 

The effect of ground conduction is described in R.ef. 6. The effect is es- 
pecially serious for,concrete runways containing steel reinforcing rods. Refer- 
ence 6 recommends using a low carrier frequency ( 165 Hz) to minimize the 
distortion effect due to ground conductivity. 


Nonlinear effects are most troublesome in the lateral deviation equation 
where the linear relationship between Ad and Az is relied upon in Eq. (76b). 
By restricting the use of the detector to rollout and turnoff, we are assured of 
operating within a narrow range of vertical displacement ( ± 2 ft. ) due to strut 
compression and pitch oscillations. Thus, we can expect to operate in a linear 
range durii^ those phases of the flight where the magnetic cable is required. 


Finally, phase modulation can be eliminated as a serious error source 
by restricting the use of the detector to small range in Az . In general, none 
of the references indicated any difficulty in obtaining a good A^ signal. 

For the purposes of simulating the error sources, the following equa- 
tions were used to generate realistic magnetic measurements of and 








V. 


NONCRITICAL FLIGHT PHASE NAVIGATIONAL AIDS 


In addition to the MLS Radar altimeter and the magnetic leader cable, 
which are required for the flight critical phases of landing approach, rollout, 
and turnoff, we include VORTAC range and bearing devices and a baro-altimeter 
as additional navigational aids for the noncritical phases of the aircraft flight 
regime. These instruments have relatively large bias errors. However, since 
they are used only for navigation between way points away from the terminal area, 
they can provide adequate control for the gyro drift and accelerometer scale 
factor error sources in conjunction with a Kalman filter. 


A. VORTAC Range and Bearing 

Let the VORTAC station coordinates be , y^ , in a level frame 

at the VORTAC station. The axes are: x^ due north; y^ due west; and 
vertical, positive upward. Let the aircraft position in the VORTAC coordinate 
system be x , y , z . In the measured VORTAC range, we have 

3i Si S, 

■•v = * ’’r 

and for the measured VORTAC relative bearing 

y ~ y 

= 'I' - = tan ^ — — — + b_ + £_ (79) 

V x^- x^ ^ >1^ 

To complete the triad, we obtain the barometer measurement of altitude 
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In the above, , b^, and b^^ are biases modeled as exponentially cor- 
related random errors and £ , £ , , and £, are white noise. The biases are 

r ijr h 

modeled as follows: 


d . r 

— b .+ a 

dt r T r y T ’r 
r r 





— b = 
dt ^ . 


^ or 


(81) 


— b = 
dt h 


\ rr 

T, ^ °h V T ^h 


where t , T, , and t, are the correlation times; a , a,, and a, are the 
r > 1 ^ h r 'I' h 

variances; and 77 ^, and r/j^ are white noise samples. (See Table 5) 

To complete the necessary equations for the Kalman filter, we derive the 
partials of the three observables with respect to the 15 -element state. For the 
range partials,' we have ‘ 
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Variable 

Time 

i , Variances 

Constant 

bias 

observ. 

Range 

900 sec. 

300 meters 

40 meters 

Bearing 

900 sec. 

0. 7° 

0.7° 

Baro-altimeter 

100 sec. 

40 meters 

40 meters 











The other 11 partials are zero. For the bearing partials, we have 


6x 

a 






5A^ 

^^a 




(82b) 


bA 'I' 





The other 12 partials are zero. 


For the barometer partials, we have 




(82c) 


The remaining 13 partials are zero. 
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VI. GUIDANCE AND CONTROL 

The Langley TCV program utilizes a modified Boeing 737 aircraft for its 
research. It contains an automated guidance and control system for landing approach 
and rollout. The simulation of the guidance laws were supplied by the Langley 
Research Center. These laws utilize ILS navigation aids. In adapting the B-737 
aircraft for the Kalman filter and the complementary filter utilizing MLS navi- 
gation aids and a magnetic cable for rollout and turnoff, certain changes had to 
be made to the Boeing Autoland gvii dance laws. The changes are outlined below: 

1. The Boeing complementary filters were removed from the 
flight critical glideslope and localizer control laws. 

2. Since no requirement for turnoff guidance existed in the original 
B-737 TCV aircraft, the Boeing rollout laws were modified 

to permit a uniform guidance law to hold for both rollout and 
turnoff. 

in general, except for the above changes, the control laws, gains, etc. used in 
the Boeing Autoland equations were maintained. 

This section contains guidance and control laws for the following phases: 

1. Flight critical lateral control for oncourse landing approach. 

2. Flight critical vertical control in landing approach for capture 
and flare. 

3. Rudder control for landing and decrab. 

4. Lateral control for noncritical level flight and coordinated 
banked turns. 
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5. Vertical control for noncritical level flight and descent. 

6. Steering laws for rollout and turnoff. 


A. Flight Critical Ijateral Control Guidajice 

The flight critical lateral guidance contains two control functions, localizer 
engage and oiicourse. These are triggered by deviation criteria computed from 
data supplied by the navigation filter and the desired glide slope path. 

The computations required for these control modes are 


■ RANGE = 'C(x_'- x)^’+ y^+ z‘ ] ^ 


V 


^ ■! y., ^ . ;. 18Q 

“ RANGE ir 


r\ 


-V-{ 


3 


range... 


z- h. 


(83) 


V = (x+y+z) 


0 =0 

, deg .. . . . 7T.. 


deg 


= P 


180 

n 


where h is the extended landing gear height (see Fig. 3). 

VJ . 


67 



The triggers for localizer engage and capture are 


L.E. 

is 

TRUE 

if 

lr?!<2.5° 


O/C 

is 

TRUE 

if 


1 1 < . 027°/sec. , 





and |0^^gl 

< 3.° 


Two block diagrams describe the flight critical lateral control mode. 

The first, Fig. 12, generates the aileron roll command signal as a function of 
the lateral deviation from the glide slope descent path. The second, Fig. 13,- 
is the flight control lateral inner loop which generates the aileron command to 
the aileron servo. 

It should be noted that the inner loop logic contains an input from the rud- 
der deerab command, . The guidance logic for this input is given in 

RDC 

Fig. 15 covering the rudder command signals. 

B. Flight Critical Vertical Control Guidance 

The flight critical vertical control guidance covers three modes:: capture,' 
capture plus 10 seconds, and flare. As in the case of the flight critical lateral 
control guidance, the inputs for these modes are generated by the navigation filter 
equations to control the vertical deviation from the glide slope trajectory. 


The equations for the required inputs are given below: 



L. E. = 25. 

L. 

O/C = 10. 

O/ 

Flare = 5. 

C2 


D1 


D3 








^LIM 

O/C = 4. °/sec 

L.E. , not O/C = 7. /sec 


C5 = 20. 

C6 = 1, 

C55 = 1 

D4 = 5. 

RD = 


775 

183 

7 


Figure 13, Flight Critical 
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Inner Loop 







( 84 ) 


y = glide slope flight path angle (-3, ) 


h 

- z 

h 

- z 

RAN = 

[(V 

^GS 

sln-1 [ 

*^deg “ 

. 180 
^ 77 


.2 .2 Jl 


{ X - x)sin y + z cos y 

Cj 

RAN 


where x is the glide slope intercept with the runway ( x = 304. 8 meters ) 
G G 


The triggers for capture, capture plus ten seconds, and flare are given 
below; ’ 

C (capture) is TRUE if 1"^ -108° 

CPIO- (capture plus 10 seconds) is TR.UE ten seconds 
after C is actuated 

F (flare) is TRUE if . 416( h^+ 15. ) + 26 h s o , 

The block diagram for the flight critical vertical control guidance is 
shown in Fig. 14. The output of this control law is the elevator command. 


C. Rudder Control for Landing and Decrab 

The rudder control guidance contains two functions; first, to provide 
yaw damping during all phases of aircraft flight, including rollout and turnoff; 
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Figure 14. Flight Critical Elevator Command. 





























secondly, to provide the rudder command for decrab immediately prior to 
rtouchdown In the event of a cross wind. 



The trigger for decrab is given by 

DC (decrab) is TRUE if O/C is true and h^ < 45. 7 meters 

The block diagram for the rudder control command for yaw damping and 
decrab is shown in Fig. 15. 

It should be noted that the rudder decrab command signal, , 

KUO ' ” 

required in the flight critical lateral control guidance loop is generated in Fig. 15. 
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Figure 15. Rudder 
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A2 = 

.3 
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• 45 
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' .34 

R3 = 

2. 22; 
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’’aw Damper and Decrab Command 









D. 


Noncritical Aileron Command Guidance 


Lateral guidance laws for trajectories between way points and for banked 
turns are governed by noncritical aileron command guidance. For banked turns, ; 
the lateral guidance requires the following information: 


- turn radius 
T 


SIGN 


- { 


+1 left turn 
-I right turn 


X , y = X and y coordinates of the center of the 
w w 

turn circle in runway inertial coordinates 


CRTE (cross track error) = SIGN 


R - /(X - 

T V ' 


v2 ■ - ,2 

X + y-y ) 
w ' w 


TANGE (track angle error) = - 


(y-y )y + (x-x )x . 

' w w' 180 


(x-x^)y -(y-y^)x 


TT 


D 


desired constant turn speed 
2 




des 


tan 


V 

V G R_ / 


180 

7T 


deg 


: -180 


: The block diagram illustrating the noncritical flight aileron command is 
shown in Fig. 16. 



C2 = (C3)V7. 12 . g 

C3 = .68 -(.00107) (.3048) [Airspeed (M/S)] ACOM 

C4 = 5. 

C5 = 1.82 

QPOT= .22 + 126.5 (.3048)/AIBSPEED(M/S) 


Figure 16. Noncritical Aileron Command. 









E. Noncritical Elevator Control Guidance 

Vertical control during the noncritical phases for trajectories between 
way points and in banked turns is illustrated in Fig, 17 for the elevator command. 

The information needed to supply the control equations for a banked turn 
along a descent trajectory for a flight path angle of -3° is 


=. altitude at the start of the turn circle 


x( t \ X X “ X 


h =0 
c 

_ .180 
Q, = q — — 
deg ^ IT 




180 

TT 


(87) 


F. Airborne Throttle Command . , 

Thrust control is designed to maintain a given airspeed and to compen- 
sate for the loss of lift during a baJihed turn. 

The variables needed to operate the throttle command guidance are 
given below: 
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f(0) 

1 

- - ( 1 - cos 0) 


V(m/s) 

airspeed in meters per second obtained from 
the airspeed indicator 


AV(m/s) = 

V - V, 
des 

(88) 

VG 

( X + y ) 


DTVG 

(ix +yy)/VG , 



The block diagram for the throttle command is shown in Fig. 18. 

G. Rollout and Turnoff Guidance ; 

Rollout and turnoff guidance encompasses rudder and/or nose wheel steering 
commands, thrust control commands for speed acceleration and taxiing, and finally, 
brake pressure commands for speed deceleration and taxiing. 

Rollout guidance is activated as soon as the main gear preload is exceeded. 

At this instant, two commands are initiated. First, the throttle is manually set 
into the idle position. This action is accomplished manually in order to prevent 
accidental actuation of the throttle into idle detent prematurely. Secondly, a two- 
second (2 sec.) slow-out and slow-in is initiated to bring the inflight, pre-touchdown 
control surface guidance commands for the aileron, elevator, and rudder to zero 
and to raise the post -touchdown rudder /nosewheel steering commands to full signal 
strength. The mixing takes place within two seconds and is illustrated in Fig. 19. 
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If AP<0. , FX = AP 

If AP > 0. , FX = I CMX - .1 - (i THRUST - KO)/K3] AP 


Figure 18. Throttle Control 



























The rollout and turnoff control is unified into a single guidance law by 
requiring the aircraft to follow a single magnetic cable. A three-coil magnetic 
pickup, mounted below the aircraft fuselage generates two signals; one measures 
the lateral deviation of the aircraft from the cable, and the other measures the 
yaw deviation of the aircraft from the cable.. The signals necessary to provide 
damping stability for the lateral and yaw deviations are generated analytically, 
in a manner similar to the banked turn logic in the lateral noncritical flight 
guidance equations. By requiring that the magnetic cable consist of a sequence 
of piecewise straight portions, connected with circular arc segments, we are 
able to specify the desired ground velocity vector and the desired, yaw rate for 
both the straight line and circular segments. By using the aircraft estimates of 
the ground velocity vector and yaw rate, we generate the deviations for use in 
guidance equations. 

In order to communicate the start of each straight line or circular seg- 
ment, a radio signal is given to the aircraft at the beginning of each segment. 

The information contains the following: 

Circular Segment ‘ 

1. radius of turn ; 

/ 

I . , ' 

2. sign of turn (+1, right turn; -1, left turn) 

3. X and y (coordinates of turn center) 

4. signal to start turn 

Straight Line Segment 

1. yaw angle of segment 

2. and y^ (coordinates of start of segment) 

3. signal to start straight line guidance 



. .The variables required to operate the rollout and turnoff guidance equations ‘ 
are listed below: ; 

For Circular Turn: ' ■- ,■ 


VG 


1 

2^ 


R 


des 


■= (x^+y-) 




Ay 


SIGN 


(x-x )x+(y-y )y 


[(i-x )^+(y-y f ] 
, w w 


1 
. a 


(89a) 


R, =rM 

deg n 


For Straight Line: 


Ay = cos y sin >1?^^ x. 


R, . = 0. 

des 


( 

i •„ 


The logic for the uniform rollout-turnoff guidance is illustrated in the block 
diagraih shown in Fig. 20. 


(89b) 


The logic for the brake pressure control is designed to first assist the 

reverse thrust in decelerating the aircraft to the coast speed (V ) , -and 

■ - i-. • f ' coast' ; 

secondly, to maintain a somewhat lower taxi speed, V . . During the first, 

• . t3.Xl 


phase. 


80% of maximum brake pressure is available to the brake control logic. 












During the second phase, only 40% of maximum brake pressure is available 
for, control. An anti-skid device maintains a modulating control to prevent 
wheel lock and skidding. 4 

The logic for the brake pressure is shown in Fig. 21. 

. The rollout thrust command is illustrated in Fig. 22. The throttle is 
held in the idle detent for two seconds. Following that, the thrust guidance 
calls for the specified reverse thrust. When the aircraft reaches a ground speed 
of V , the throttle control is once again brought to the idle position and’ ' 

COaiSu 

remains there throughout the taxi phase. 

. ■ -j . • L 

,'4 ■ . , 1 . 

- ' . ' *' 5 *' ■ , \i 

H'. Requirements for GE Whole-Word fcomputer to Execute Navigation- and 

. ,1 • ' i'. 

Kalman Filtering During Rollout - ' , 

. . This section contains estimates of the timing.and storage requirements 

for the General Electric MGP-701 computer to execute the navigation and Kalman 
filtering tasks during the touchdown and rollout sequences. Since the navigation . 
and filtering functions are, programmable as separate and dis'ti’hct'.mbdules; esti- 
mates of their timing and^torage requirements may, be made separately, and'- 
the resulting estimates may be summed to, obtain total estimates for timing and , 
storage. < .1 ^ 

• Different methods were used in obtaining the estimates for the navigation, 
function and for, the Kalman filtering function. 

• FORTRAN code of the rollout navigation algorithm was used in making the 
storage and timing estimates. FORTRAN operations were translated into; equivalent 
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MCP-701 operations in estimating both timing and the program-instructions 
component of storage. 

A more direct approach was used for the Kalman filter. The Kalman filter 
is already programmed and running in single precision on the. Sperry 1819A com- 
puter at the Ames Research Center. Therefore, the timing estimate is the actual 
1819A running time scaled to reflect the relatiye performance of the MCP-701. 

The storage requirement for the filter is essentially identical to that of the 1819A. 

Kalman Filter Requirements ; 

i - • . • 

t " ■ 

To give the estimates of the filter storage and timing requirenients some 
meaning, a brief description of the Kalman filter is ih order. It is a three-axis 
Kalman filter with segregated horizontal and vertical channels, as depicted in 
pig. 23'(for the Sperry 1819A computer). For the application under consideration, 
the MODILS measurements would be replaced by the equivalent MLS measurements, 
and the TACAN and MLS measurements would probably not be used simult^eously.' 

Attitude and heading data and three body-mounted accelerometers are used 
to give the acceleration in runway reference. Existing 1819A software performs 
these calculations at a 20 Hz frequency. Biases in the accelerations (estimated 
by the x-y and z filters) are added to the accelerations and then integrated in 
double precision to give the three-axis velocity vector, which is in turn integrated 
in double precision to give the three-axis position vector. This is the only double- 
precision logic and is performed in the block Navigation Equations in Fig. 23. The 
estimated state, vector coming from this block includes the following: 

3 -component position vector with respect to runway 
3-component velocity vector with respect to runway 
3-component acceleration bias vector 




Figure 23. Block Diagram of Overall Filter Mechanization 














2-component estimate of winds (:x and y ) 

1-component estimate of TACAN range bias 
1 -component estimate of TACAN bearing bias 
1-component estimate of baro-altimeter bias 

The position and velocity components of the estimated state are updated from the 
acceleration data at the 20 Hz frequency. All components of the state are up- 
dated from the x-y and z filters at .'667 Hz , This low-frequericy update causes 
discrete changes in the estimated state, which are likely to be objectionable if this 
estimated state were used for guidance display and control purposes. The block 
labeled Smoothing Logic in Fig. 23 filters the discrete changes from the estirhated 
state so that the smooth state estimate output will be free of objectionable jumps. 

The x-y filter accepts the true airspeed, TACAN range and bearing, and 
MODILS range and azimuth measurements. This filter has 10 states, as follows: 

X - position along the runway 

y - position normal to the runway , . 

X - velocity along the runway. 

y ■ velocity normal to the, runway , , 

- along runway component of acceleration bias 
.. . , a^ - normal to runway component of acceleration bias , , 

b - bias in TACAN range measurement 
, ~ ^ . TACAN bearing rneasuremept 

Wx .wind along the runway . .... 

w^ - wind normal to the runway , . 

The z filter accepts altitude measurements from the barometric altimeter 
and radio altimeter and altitude measurements calculated from the MODULS eleva- 
tion data. This filter has four states: 



z - vertical component of position 

z - vertical component of velocity 

a - vertical component of acceleration bias . 
z 

b, - bias in the barometric altitude 
h 

All raw measurements are accepted at a 10 Hz frequency. Preprocessing 
algorithms perform residual and partial-derivative calculations as well as validity 
. checks, and the resultant good data is accumulated for a period of one-and-one- 
half seconds. The accumulated information is then transferred to locations for 
processing by the square-root Kalman filter algorithm. Once the filter has esti- 
rnated the error state for this information, the logic updates the estimated state 
and performs the smoothing. This concludes the brief description of the Kalman 
filter. 

The total storage (instructions and common variables) of the filter, which 
is mostly in single precision, is 2,825 words on the Sperry 181 9A computer. If 
the filter were entirely in double precision, it would require 3, 053 words. These 
storage requirements are essentially valid for any mini-computer, since the 
assembly-language instructions for various mini-computers are very similar. 

The filter consumes 18.3% of 1819A running time. The MCP-701 is felt 

to be at least twice as fast as the 1819A. In particular, double -precis ion additions 

' 

are approximately twice as fast and single-precision multiplications are about four 
times as fast. Therefore, it is estimated that the filter will consume about 10% 
or less of MCP-701 running time if implemented mostly in single precision. A 
consideration of the operations necessary to implement a double-precision multiply 
function on the MCP-701 has yielded the conclusion that a worst-case factor of 4 
would result if the filter were programmed entirely in double precision. It is esti- 
mated, therefore, that the filter would consume about 40% or less of MCP-701 
running time if implemented entirely in double precision. 
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Navigation Requirements 


' 'A detailed investigation of each phase of the rollout navigation algorithm 
was performed to determine the timing and storage requirements. FORTRAN 
operations were translated into equivalent MCP-701 operations (or sets of opera- 
tions). Timing estimates for each phase of the rollout were then estimated by 
directly counting the required equivalent MCP-701 operations and making use of 
the information in Table 6. 


The rollout sequence and timing, estimates are presented in Fig. 24. The 
times are given in microseconds and include a 50% "safety factor" over the values 
actually Computed (i. e. , a multiplicative factor of 1. 5). Timing estimates are 
given for the assumptions of all-single-precision arid all-double-precision coding, 
and percentages of real-time consumed are given in parentheses assuming a 20 Hz 
frequency. Two situations may occur during rollout: 


; : a) braking is completed before the turn begins;,' * , . 

b) braking continues into the turn 

The latter easel results in theTargest percentage tiJriC load, of the -rollout navigation 
algorithm on the computer ( 2, 5% in single precision and 8. 6% in double precision). 


, The rollout navigation algorithm is estimated to require 702 MCP-701 in- 
structions. The size of COMMON (where variables and constants are stored) is 
estimated to be 146 words if single-precision coding is sufficient, , and 277 words 
if double precision is required. These values include a 50% "safety factor" over 
the values actually computed. The total storage estimate for rollout navigation is 
obtained by summing the instructions and COMMON values. In the, single-precision 
mode, a total of 848 words are required, and in double precision 979 words are 
required. 


92 



Memory Addressable 
Instruction 


MEMORY ADORESRABLE INSTRUCTIOM FORMAT 



MODE VARIATION 

TK« mofh v«ri«tiora for mamory addrvaaOt* intowcttons 
(Mint tha machod for Qtnarannf tha affactiaa aOtfras 
location of diu within mamory. 


MOOC 

E^^ECTIVE AOOftESS 
GENERATfON 

i^TlME 
(m «ae) 

m 

n 

a 



iD 

o 

o 

Currant Paja 

0 

a 

o 

• 

2aro Paga 

0 

o 

a 

o 

Indiract 

1 ■ 

n 

' 

' 

2aro Paga/lndirect 

1 

n 

o 

a 

Indai 

a2. 

n 

o 

' 

- 2^0 Paga/lnda« 

&2 

1 

1 

O 

Irtdiract/IndaK 

1 

1 

t 

s 

2aro Paga/lndiract/lrtdaa 

.1 


LOAD INSTRUCTIONS assemile* EORUAT 1 


UNE 

CODE 

INSTRUCTtON NAME 

TIME IM MCj 

LCU 

OB 

Load Uppar Ragittar 

2.2 

LOL 

10 

Load Lowa* Ragittar 

2.2 

LOX 

• 18 

Load Inbai Raginar ' 

2.2 

LOO 

20 

Load Ooubta 

X2 

LPK 

28 

Loai Packad 

2.2 


STORE INSTRUCTIONS 


STU 

30 

Stort UpCMT Ragictar 

2.2 

STL 

38 

Siort Lowar Ragittar 

2.2 

STX 

40 

Slora Indaa Ragittar 

2.2 

STD 

48 

Stora Ooubta 

3.2 

SPK 

50 

Stora Packad 

2.4 


ARITHMETIC INSTRUCTIONS 


AOU 

sa 

Add to Uppar Raginar 

2.2 

ADM 

60 

Add to Mamory 

X2 

ADO 

68 

Add Ooubta 

X2 

SBU 

70 

Subtract from Uppar Ragiatar 

2.2 

S8D 

78 

Subtract Ooubit 

X2 

Mpy 

60 

Multiply. 

5.4 

OlV 

88 

Oiwida 

6.8 


LOGIC INSTRUCTIONS 


ANU 

90 

And (0 Upper Ragittar 


ORU 

98 

Or to Upper Hegiftar 


EXU 

AO 

Ezdusivt Or u> Uppar Ragittar 


CPU 

AB 

Compare to Upper Ragittar ' 

1 


BRANCH INSTRUCTIONS 


JMP 


Jump Uncckodiiioiutlv 

1.4 

JMS 


Jump to Subroutine 

2.7 

* MOS 

CO 

Mamory Dacramant and Skip 

3.2 


Non-Memory Addressable 
Instruction 


NON-MEMORY ADDRESSABLE INSTRUCTION FORMAT 


IS 14 13 12 11 

MOOC 
18 • 8 

F2 

7 6 8 4 

n : 

3 2 10 

nTTT 

M 

_L 

Ml 

INSTnuCTION 



REGISTER COOES 


MNE 

CODE 

REGISTER 

UR 

1 

Uppar Ragittar 

LR 

2 

lowv Ragittar 

MR 

3 

Indaa Ragittar 

SR 

4 • 

Stitua Ragittar 

UL 

C 

Uppar aryf Lower Ragitt«' 


INDICATOR CODES 



1 CODE 1 

indicator 

SA 


Sign of Adder 

ZA 


Zero Adder 

- AF 


Arithmetic Fault 

CF 


Computar Fault Oitablf . 

DR 


Device Ready 

10 

^B^H 

Input/Output Ortabla 

LA 


Logic Iridicatar A 

L8 

mm 

Logic Irtdicator 9 


BRANCH INSTRUCTIONS 


MNL 

COOL 

(HyTMUCtIUN NAMl 

f7 

n 

IiUC 
Im iki 

RA>>N 

ai 

AaW ft U «> ur> 




HSJrr 

CA 

a Skiu on 


C«n»i4ni 

1 ft 

RAtf 

CC 

AtMl a Skai «•» 


Ci^urt 
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RSiiP 

C£ 

SwDtraei ft Skip on Rm 

R«V>»iar 

CoMBUnt 

l.« 

SPI 

bd 

^■0 f erw«ro on Irtf . 

lne««nr 

Cenauni 

1.4 

sai 

01 

Skip ft«*>»areon »gd 

lne««IOT 
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1.4 


07 
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03 
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VH 

04 
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S9M 

Oft 

Skb iKhwaiewt Me li^ 
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1.4 


SET/ 

^CLE 

AR INSTRUCTIONS 



pi 


Indiuia* 

InoMWr 

Net Uwi 
Noi Uad 

1.4 

1.4 

REGISTER INSTRUCTIONS , 

RTF 

2R0 

IMV 

AIS 

SO 

El 

E2 

E) 

E4 

Roaunr Trmdar 
Zve RByrnar 
Caiixiomiru RoyiMr 
irwoo RopMiar - 
AtMOiwIo V«hM el Rogitw 

Frem Rayitur 
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Fiom RtgittB* 
From RB9>tttr 
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re R«giiar 
re Raetaar 
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\M 
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tA 
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Eft 
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EA 

R«ai Seacae Dwiea 

Not Uad 
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EC 

TtAa Oau le Uoew Reyaar 

Net Uad 

tA 
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EO 
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tA 
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Ei 
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tA 

»L 

EF 
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SHIFT INSTRUCTIONS 


I 

FO 

FI 

F2 

FJ 

FS 

FA 
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4 
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4 
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Table 6. GE MCP-701 Computer Instruction Execution Times. 
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Figure 24. Rollout Navigation Timing Estimates for the MCP-701 Computer 



Total Requirements 


The total Kalman filter and rollout navigation requirements (with a 50% 
safety factor included in the navigation component) are obtained by summing the 
separate requirements discussed above, assuming the worst case for the rollout 
timing estimate. In single precision, the total requirements are 3,673 words 
of storage and 12. 5% of real time; in double precision, the total requirements 
are 4, 032 woirds of storage and 48. 6% of real time. 
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APPENDIX A 


B-737 

Aerodynamic Force and Moment Gdefficients 

and 

Dynamic Response 
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I. 


AERODYNAMIC COEFFICIENTS 


The aerodynamic coefficients are functions of the velocity of the aircraft 
relative to the air stream, control surface deflections, altitude and altitude rates 
of the aircraft. 


We define the following quantities: 


c 

b 

S 


P 

6 . 


STAB 

^SP 

^SPSUM 


q: 


mean aerodynamic chord 
wing span 

wing reference area 
air density 
elevator deflection 
aileron deflection 
stabilizer deflection 
rudder deflection 
aileron spoiler deflection 
symmetric spoiler deflection 
fuselage offset angle 


The velocity vector of the aircraft relative to the airstream in the runway 
coordinate system is given by 


V = R - {W+G + S } 


(A-1) 


The relative velocity vector in the body system is given by 

V, 


3B 


IB 
''2B 


BI 


{v } 


(A-2) 
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The airspeed is the magnitude of the vector V . 


(A-3) 


The dynamic pressure is given by 


1 2 
q = 2 pv 


(A-4) 


The angle of attack is given by 


a = sin 


3B 


/v ^+v 2 
V IB 3B 


or 


O' = tan 


-1 '^SB 


(A- 5) 


IB 


The aerodynamic coefficients are functions of the angle of attack relative 
to the mean chord line 


a = a + a 


(A-6) 


The side force angle of attack is given by 


^ = 


2B 


(A-7) 


The effective drag coefficient is given by 


C = C +C a+(C +C a)f> +C 

D 0102^^ ' D3 D4 '^SPSUM D6 


+ + CdioO!)6sp+ Cd 


(A-8) 


GEAR 



The effective lift coefficient is given by . 


^L1 ^L3 2v^“ ‘^gust^ ^L4 2v ^ ^ \ust^ 


'*' ^L5 ^STAB ^ ^L6 ^ ^L8 ^SPSUM 




GEAR 


The effective side force coefficient is given by 


Cy - Cy^ + Cyg ^ ( P ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 


■*■ ^vo 7 T"(r -psina + r cos a)' 
Y3 2v gust 


(A-10) 


+ (Cy 4 + Cyg o;)6gp + Cyg 


The aerodynamic roll moment coefficient is given by 


B + — (p :, + p cos a + r sin a) v 

SI S2 2v ' gust 
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The pitch moment coefficient is given by 


^M1 ^ ^M2 ^ ^ ^M3 2v ^ ^ ^gust ^ ^ ^M4 2 v ^ ^ ^gust ^ 


^M6 ^STAB ^M7 ^MIO 1 ^ I ^ 


Mil 


M = qSc ^ 
o 


^^M12'^ ^MIS^^^PSUM"^ ^M14 I 


■^C Mg-.25)C 

GEAR 


In order to compute the yaw moment coefficient, we require the lateral 
acceleration in the body reference frame. We have 


ii 

V 

w 


> = L 


BI 


{M 




The yaw moment coefficient is given by 


Ng = qSb < 


b V b 

C-T 1 0+ C-.„ + C.,„ t;- ( P . + P cos a + r sin a ) 

N1 N2 2v V N3 2v gust > ; 


+ .-psina+r cos a ) + C.,„ 5„_ 

N4 2v ' gust. ‘ N6 '^SP 


+ (C +C a)6 +C 6 +-(g*-.25)C 
' N7 N8 ' a N9^R r ’ Y 

0 


The numerical values of the coefficients for the 737 aircraft were fur- 
nished by NASA Langley Research Center and are listed below: 


> (A -12) 


-13) 


(A-14) 



AIRCRAFT AERODYNAMIC COE FFICIENTS 


Drag Coefficients 


'"di 

- 

.185 

^D2 

S 

.9225 

^D3 


.0329 

^D4 


-.1067 

^D5 


0. 

^D6 

= 

.0075 

^D7 


.0573 

00 

Q 

O 

— 

.019 


=r 

.043 

^DIO 


-.308 


Side Force Coefficients 



-1. 564 

^Y2 

. .4871 


.189 

O 

II 

-.0487 

. II 
O 

-.1067 

°Y6 = 

.4383 


Lift Coefficients 



= 

1.36 

^L2 


6.9328 

^L3 


-8.0 


= 

8.208 

^L5 


.963 

^L6 

- 

.464 

'^L7 

' = 

-.4584 

^L8 


-.0258 

^L9 

r; 

-.058 

^LIO 

= 

93. 

''m 

= 

.008 

^L12 

= 

-.625 

^L13 


.7735 

Roll Moment Coefficients 

Si 


-.3152 

S2 

= 

-.645 

S3 

= 

.395 

Si 

— 

.0967 

Ss 

zz 

. 0817 

Se 

--- 

.0816 

St 


.06 

Ss 

= 

-.1899 




c 

Me 

c 

M 7 

C, 

Ms 

C 

Ms 

c 

Mlo 



M12 


C 

M13 

C 

M24 


437 

6545 

'loss 

•4842 

"• 34 


•OO3 
'•02 
'• Ssg 
•Ms 




•(^215 

•Oil 


Ns 

c 

Ns 


• 0343 


> 


The geometric parameters for the 737 aircraft are listed below: 


c = 
b = 
S = 


Mass = 
^XX "" 

lyy " 

^ZZ " 
^XZ ~ 


3, 41376 meters 
28. 3464 meters 

2 

91. 045 meters 
1. 524 meters 
36287. 5 kg 

2 

142572 newton meters 

2 

328148 newton meters 

2 

490738 newton meters 

2 

21351 newton meters 


(11.2 ft.) 

(93 ft.) . 

(980 ft.^) 

(5 ft.) 

(80,000. lb.) 

(^5,000. lb. ft.^) 
(794,062,5 Ib.ft.^) 
(1487,500. lb. ft.^) 
( 51,667. Ib.ft.^) 



^ / 
/ 

II. DYNAMIC STABILITY COMPARISON ' 

In order to check the accuracy of the dynamic model of the aircraft, the * 

response of the plane to an elevator impulse ^d a rudder impulse was tested. 

The response curves are, shown in Figs. A1 and A2. These r^hs were made with." 
the automatic control equations inoperative in order to obtain the period, damping, 
and the decay time to one-half amplitude of the plane for the short period, phu- , 
goid, and dutch roll mptions, at an airspeed of 66. 88 m/sec. (130. knots) and : 
an altitude of 426. 7 meters (1400 ft. ) 

A comparison of these results was made with 737 aircraft characteristics 
supplied by Langley Research Center. The results, shown in Table 7, are in ' 
good agreement. 
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Dynamic Response to Elevator Impulse 


TIME 



^ 4 e e 10 12 14 18 18 20 X 10* 



4 6 B 10 12 14 18 18 20 X 10‘ 



Figure Al. Phugoid and Short Period Modes 


SDESI.rP ; YAWRATE RUDDER YRf ROLL, 



TABLE 7 


DYNAMIC RESPONSE COMPARISON — B-737 


Simulation Model 


Mode ■ 

. Period 

,, Damping 

: ’^1/2 

(sec. ) 

Coefficient 

(sec.) 

Phugoid 

39.0 

. 098 

43, 7 

Short Period 

5.0 

.25 

4.84 

Dutch Roll 

5.0 

.055 

9.99 


B-737 Data 


Phugoid 

34.0 

.084 


43.7 

Short Period 

5.39 

. 41 


5,39 

Dutch Roll 

■ 4.81 

" : ' ' ^057 ' 

., . . .. . . 

" 9:31 



APPENDIX B 


Computer Simulation Results for 
Rollout and Turnoff 



I. 


DESCRIPTION OF TEST DATA 


This section contains the results of computer runs of the ALERT program 
in the form of plots. Although the B-737 aircraft simulation program is capable 
of executing automated banked turns and landing approaches including capture, 
oncourse, flare, and decrab, the main purpose of this study was to develop and 
test a guidance law for rollout and turnoff. The computer plots contained in this 
section illustrate the terminal portion of each landing from touchdown to the end 
of turnoff. Turnoff, in the context of this study, is defined to be completed when 
the aircraft is at a fixed desired taxi speed, at a fixed constant heading ( 30° yaw 
relative to the landing runway centerline), and at a distance of 91.44 meters 
(300 ft. ) from the runway centerline. 

In order to reduce the computer run time to complete the study for a 
range of wind conditions, runway conditions, and taxi speeds, the entire air- 
borne program was run for different wind conditions. The aircraft variables 
at touchdown were stored, and a set of rollout and turnoff runs were made 
at different taxi speeds and runway conditions (dry and wet) for each fixed landing 
condition. 

A description of the plot variables and units is given. Each landing case 
contains eleven plots. The horizontal variable is always time from touchdown 
in seconds. The first plot is the aircraft thrust in newtons. (One pound of force 
is equal to 4. 4482216 newtons.) The specified reverse thrust for this aircraft is 
62275. newtons (14,000.#) The steady state idle thrust is 6672. newtons (1500.#) 
The thrust is manually set at idle and, after a two-second slowout, the thrust 
logic calls for maximum retro thrust which continues until the aircraft reaches the 
required V speed, at which time the thrust logic calls for idle thrust through- 

COaS L 

out the taxi period. 



The second plot is the automatic braking pressure plotted in percent of 
maximum allowed pressure. 100% of maximum allowed braking pressure 
corresponds to a friction coefficient, of .4. Thus, the automatic braking 

logic is not designed to command the full capabilities of the tire on a good dry 
runway. Zero brake pressure corresponds to the condition where only rolling 
friction is operative. Maximum brake pressure is available during deceleration 
to V . Following that, only medium brake pressure (fj, = .2) is available 
to the automatic braking logic. During landings on wet runways, the anti-skid 
device moderates the pressure to that allowable to forestall wheel lock. 

The third plot is the forward acceleration in the body system, u . The 
2 2 

units are in meters/sec. . l.Og is 9. 817 m/sec. . For rapid turnoffs, a 

2 

maximum deceleration of 5 m/sec. is considered within the limits of passenger 
comfort. 

The fourth and sixth plots are the x and y components of the winds plus 
shear and gust. These are measured in meters/sec. 

The fifth plot is the aircraft ground speed. This is measured in meters/sec. 

The flat portion after the linear deceleration illustrates how well the aircraft keeps 

the desired V. . speed, 
taxi 

The seventh plot is the first plot on the second page of each case. The 
variable CRTE is the cross track error, or the lateral deviation from the 
buried magnetic cable. It is measured in meters. 

The eighth plot is the lateral acceleration in body system. It is measured 

in meters/sec. Passenger comfort requires that this variable not exceed an ab- 

2 2 

solute value of .ISg’s or 1.472 m/sec. (4.83 ft. /sec. ). This limits the maxi- 
mum taxi speed if fast turnoff times are required. 
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The ninth plot is the aircraft yaw angle. 

The tenth plot is the nose wheel steering angle. 

The final plot is the difference between the magnetic cable yaw angle and 
the aircraft yaw angle. 

For the purposes of this simulation, the magnetic cable is laid out in three 
segments; first, a straight segment along the runway centerline to the turnoff 
position; secondly, a turn with a fixed input radius making a change in yaw of 
30° heading; finally, a straight line segment at a fixed yaw heading of 30^ , 

Following the plot of the eleventh variable is a short data set giving the 
wind condition, the V and V . inputs, and the runway condition, 

OO^LSl ' 
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II. 


SIMULATION RESULTS AND DISCUSSION 


Figures Bl, B2, B3, and B4 are acceptable fast turnoffs bn good dry 
runways. Figure B5 is a repeat of Run_B3_.for a wet runway. For this case, 
because of the inability to apply braking at high speeds, the desirable taxi speed 
of 50 knots was not achieved and high lateral acceleration resulted. In order 
to produce acceptable turns at high speeds on wet runways, it was found neces-r 
sary to lower the V speed (the speed at which the retro thrust is switched 

C03.St ■ r 

to idle) to exceed the desired V . speed by 5 knots. With this change, ac- 
ceptable fast turns are achieved on wet'-rvinways. , Figure. B6 is an illustration ! 
of a good 6. 0 knot turn on a wet runway. 

Figures B7 through B14 are all acceptable high-speed turns on wet run- 
ways for different wind conditions. ‘ . 

Figure B15 is an illustration of high-speed turnoff on a good dry runway, 
replacing the Kalman filter with the complementary filter in the navigation loop. 
Examination of the- cross track error in Fig. BI5 shows-that the complementary 
filter does a good job until the turn is coinpleted, at which time it ends up with . 
an error in lateral deviation of 5 meters as compared with the Kalman filter 
deviation of only 2 meters . 
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Figure Bla. High-Speed Turntoff on Dry Runway 
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hind Y.n/5: VG.tt/S , WIND X.n/5 UpOT..n/S/5 FERCNT SRK THRUST. N 








DELYflW.OEG ^ NOSMHL.OE(j YflU.OEG VOOT.n/S/S CRTE.n 



30° Tail Wind, 5. 144 m/sec (10 knots) 
V .= 30.86 m/sec (60 knots) 

COaSt 


V = 20.57 m/sec (40 knots) 
taxi 

R^=(2500 feet) 


116 

c 


WIND V.n/s VK n/fi KINO X,n/5 UqOT.rVS/S , rWCNT BKK THRUST.N 



Figure B3a. High-Speed Turnoff on Dry Runway 
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Figure B3b. High-Speed Turnoff on Dry Runway 

30° Tail Wind, 5. 144 m/sec (10 knots) 

V 36.01 m/sec (70 knots) 
coast 

V = 25.72 m/sec (50 knots) 
taxi 

R.p= (3500 feet) 
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Figure B4a. High-Speed Turnoff of Dry Runway 
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Figure B4b. High-Speed Turnoff on Dry Runway . 

30° Tail Wind, 5 . 144 m/sec (10 knots) 

V =41.16 m/sec (80 knots) 
coast 

y . = 30.87 m/sec (60 knots) 
taxi 

R^ = (3500 feet) 
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Figure B5a. High-Speed Turnoff on Wet Runway 
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5, 144 m/sec 
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Figure B6a, High-Speed Turnoff on Wet Runway 

/ 
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Figure B7a. High-Speed Turnoff on Wet Runway 
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Figure B8a. High-Speed Turnoff on Wet Runway 
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Figure B8b. High-Speed Turnoff on Wet Runway 


150° Head Wind, 5, 144 m/sec (10 knots) 

V =28.29 m/sec (55 knots) 
coast 

V = 25.72 m/sec (50 knots) 
taxi 

R_ = (3000 feet) 
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Figure B9a. High-Speed Turnoff on Wet Runway 
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Figure B9b, High-Speed Turnoff on Wet Runway 

150° Head Wind, 5, 144 m/sec (10 knots) 

V =33.44 m/sec (65 knots) 
coast 

V = 30.87 m/sec (60 knots) 
taxi 

R^= (3000 feet) 
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Figure Bllb. High-Speed Turnoff on Wet Runway 

150° Head Wind, 7.716 m/sec (15 knots) 
V = 28.29 m/sec (55 knots) 

CO aSX 

y .=25.72 ra/sec (50 knots) 

taxi 
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Figure B13a. High-Speed Turnoff on Wet Runway 
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Figure B15a. High-Speed Turnoff oh Dry.Runway 
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Figure B15b. High-Speed Turnoff on Dry Runway 

150° Head Wind, 5. 144 m/sec (10 knots) 

V = 30.87 m/sec (60 knots) 
coast 

V =20.53 m/sec (40 knots) 
taxi 

R^=(2000 feet) 

Complementary filter used 
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